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Abstract. - We present a fully electronic analogue of coherent population trapping in quan- 
tum optics, based on destructive interference of single-electron tunneling between three quan- 
tum dots. A large bias voltage plays the role of the laser illumination. The trapped state is a 
coherent superposition of the electronic charge in two of these quantum dots, so it is destabilized 
as a result of decoherence by coupling to external charges. The resulting current / through the 
device depends on the ratio of the decoherence rate and the tunneling rates. For P^, — > 
one has simply / = eP^. With increasing P^, the current peaks at the inverse trapping time. 
The direct relation between I and can serve as a means of measuring the coherence time of 
a charge qubit in a transport experiment. 



Coherent population trapping is a quantum optical phenomenon in which the laser illu- 
mination of an atom drives an atomic electron into a coherent superposition of orbital states 
and traps it there [1-3]. Such superpositions can be "dark", in that they are further decou- 
pled from the optical fields. Brandes and Renzoni have shown how such states can also be 
formed in artificial atoms (quantum dots) through the use of laser illumination [4,5]. In this 
paper we present an all-electronic analogue, i.e. without laser illumination, of coherent pop- 
ulation trapping in quantum dots. (For an analogy in superconducting Josephson junctions, 
see Ref. [6,7]; for an analogy in single benzene molecules, see Ref. [8].) We illustrate this effect 
by considering a system of three tunnel- coupled quantum dots and show that, under proper 
bias and resonance conditions, an electron can become trapped in a coherent superposition of 
states in different dots. This state is "dark" in the sense that, due to the Coulomb blockade, 
no further electrons can pass through the dots and current flow is blocked in the absence of 
decoherence. 

The trapping effect provides a novel mechanism for current rectification, since the block- 
ing is effective for one sign of the bias voltage only. This quantum mechanical mechanism is 
distinct from mechanisms discussed previously. In particular, the classical rectification mech- 
anism of Stopa and collaborators [9, 10] traps the electron in a single quantum dot, rather 
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Fig. 1 - Schematic of the three-dot trap. The solid arrows indicate reversible transitions, described by 
the Hamiltonian fT). The dashed arrows indicate irreversible transitions, described by the quantum 
jump operator Destructive interference of the two reversible transitions traps an electron in a 
coherent superposition |<J>-) = 2~ 1 ^ 2 (\A) — \B)) of the states on dots A and B. The trapped state 
has vanishing amplitude on dot C, so that it can not decay into the right reservoir. No trapping is 
possible if the bias is inverted, because the trapped state would then decay into the left reservoirs. 

than in a coherent superposition of spatially separated states. Experiments by Ono and col- 
laborators [11] on rectification in double quantum dots likewise trap an electron in a single 
dot. The three-dot configuration requires no Aharonov-Bohm phase to trap an electron, in 
contrast to the two-dot configuration of Marquardt and Bruder [12]. Because of the phase 
coherent origin of the effect discussed here, the current that leaks through the device when 
it is blocked provides a method by which one can determine the coherence time of a charge 
qubit. 

The three-dot trap is shown schematically in Fig.^ Three quantum dots and three electron 
reservoirs are connected by reversible or by irreversible transitions. Reversible transitions 
between the quantum dots are described by the tunnel Hamiltonian 

H = T\C)(A\ + T\C)(B\ + H.c. = 2 1 / 2 T|C)($+| + H.c. (1) 

We have defined the states 

|*±) = 2- 1 / 2 (|^) ± (2) 

We consider the case that the energies of the single-particle levels \A), \B), \C) in the three 
dots are all the same (set at zero), so that inelastic transitions between these levels do not 
play a role. To minimize the number of free parameters, all tunnel rates are put equal to T . 
(The more general case of unequal tunnel rates will be considered at the end of the paper.) 
We assume time-reversal symmetry, hence T is a real number. (Since results do not depend 
on the sign of T, we will take T positive for ease of notation.) We furthermore assume that 
the electrostatic charging energy of the combined three-dot system is sufficiently large that 
the total number of electrons does not exceed one. Many-electron states are projected out 
and hence we may ignore spin. 

For a bias voltage |V| 3> T/e, and at zero temperature, the transitions from the source 
reservoirs into dots A and B and from dot C into the drain reservoir are irreversible. (Because 
of this restriction, the rectification provided by our device does not apply to the range |V| < 
T/e around zero bias.) The tunnel rates between dots and reservoirs are all set equal to T. 
The quantum jump operators are 



L a = Vt\A)(0\, l b = Vt\b)(o\, L c = Vt\0){C\, 



(3) 
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where |0) is the state with all three dots empty. 

We study the dynamics of this device by means of the master equation approach to single- 
electron tunneling [13-15], which describes not only the populations of the dot levels, but also 
accounts for quantum coherences between them. The master equation gives the time evolution 
of the three-dot density matrix p(t) in the Lindblad form [16] 



X=A,B,C 



7~1 [LxP^x lL x L x p — \pL x L x 



(4) 



(We have set h = 1.) As initial condition we take p(0) = 
We use as a basis for the density matrix the four states 



| ei > = |<& + >, |e 2 ) = |$_>, \e 3 ) = \C), |e 4 ) = |0). 



(5) 



This four-dimensional space may be reduced to a three-dimensional subspace by noting that 
the master equation Q couples only P44 and /?« with i, j < 3. The matrix elements pij with 
i = 4, j 7^ 4 or j = 4, i ^ 4 remain zero. We may therefore seek a solution of the form 



p(t) = p(t) + [1 - Trp(i)] |0)(0|, 



(G) 



where p is restricted to the three-dimensional subspace spanned by the states | e^) with i < 3. 
The evolution equation for p reads 



dp 



dt 

M = 



= Mp + pM^ + Q, p(0) = 0, 



2V 2 iT\ 

2 1 / 2 iT T/2 

1 o^ 

Q = T(l- Trp) (010 

0; 



(7) 
(8) 

(9) 



All off-diagonal elements of p vanish, except the purely imaginary P13 = —p^i- Four real 
independent variables remain, which we collect in a vector v = (Tr p, pn, P33, Impi3) satisfying 



dv 

~dl ~ 

X = 



X ■ (v - Woo), v(0)=0, 
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(11) 



The solution is 



e xt v Q 



v(t) = v 

All four eigenvalues A„ of X have a negative real part, so v(t) 



for * 



(12) 

00 and hence 



lim p(t) = |$_)($- 

t— >-oo 



(13) 
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Fig. 2 - Solid curve: Dependence of the trapping rate a on the tunnel rate V between quantum dots 
and reservoirs. Both rates are normalized by the tunnel rate T between the quantum dots. The small 
and large-r limits are given by Eq. 1141 . Dashed curve: Maximal steady state current 7 ma x in the 
presence of decoherence, according to Eq. 1191 . The two quantities / max and ea differ by less than a 
factor of two over the whole range of tunnel rates. 



This is the trapped state: it does not decay because it is an eigenstate of H. For large times 
\v(t) — fool oc e~ at , with trapping rate a — min(|ReAi|, |ReA2|, |ReA3|, |Re A4 1 ) . The full 
expression for a is lengthy, but the two asymptotic limits have a compact form, 

f 4T 2 /r if r«r, 

a { i(5- Vr7)r « o.22r ifr<r. [ ' 

If the coupling of the quantum dots to the reservoirs is weaker than between themselves, then 
the trapping time is of order 1/r. One might have guessed the trapping time to be of order 
1/T in the opposite regime T -C T, but this guess underestimates the correct answer, which is 
larger by a factor T/T. The fact that a — > when r — > 00 can be understood as a decoherence 
of the inter-dot dynamics induced by a strong coupling to the electron reservoirs. 

The full dependence of a on V and T is shown in Fig. [5J If T is increased at constant T, 
the trapping rate has a maximum of a max = 0.58 T at V — 4.35 T. 

The trapping effect requires the coherent superposition of spatially separated electronic 
states in quantum dots A and B. Such a charge qubit is sensitive to decoherence by coupling 
to other charges in the environment, which effectively project the qubit on one of the three 
localized states \A), \B), \C). We model this decoherence by including into the master equation 
the three quantum jump operators 

L^ x =T l J 2 \X){X\, X = A,B,C. (15) 

The decoherence rate parameterizes the strength of the charge noise and is taken to be dot 
independent. For a microscopic foundation of the charge noise model we refer to Ref. [17]. 
We also note that charge noise causes phase as well as energy relaxation^ 1 ) 
The master equation reads 

^ = -i[H,p}+ { L xPL x -\L\L x p~\pL^ x L x \ (16) 

X=A,B,C,<t> A ,tt>B,<t>c 

( 1 )To calculate the energy relaxation, we decouple the three quantum dots from the electron reservoirs (setting 
T = 0) and calculate dE/dt = (d/dt)Ti pH from Eq. O^J. One finds dE/dt = -T^E, so the energy of the 
three-dot system relaxes to zero with rate T^. 
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The steady-state current, 

7= lim eT(C\p{t)\C), (17) 

t — >oc 

is obtained by solving Eq. l|16f) with the left-hand-side set to zero. We find 

AeTT 2 



I 



r 2 + i4T 2 + 2rr^(i + 2T 2 /r 2 ) 

eT^ if < T, T, 

2eT 2 /r if r »r,T. 



(18) 



As illustrated in Fig. [3 the current vanishes both in the limit — > 0, because of the trapping 
effect, and in the limit — > oo, because of the quantum Zeno effect [18, 19]. The maximal 
current is reached at = 2 1 / 2 T and is equal to 

4err 2 



(19) 



T 2 + 14T 2 + 4V2rT 

4eT 2 /r if T < r, 

|er w 0.29 er if r«T. 

Comparison of Eqs. 1)14(1 and ((19(1 shows that the maximal current /max in the presence of 
decoherence is set by the trapping rate a in the absence of decoherence. For TcT one has 
exactly J max — ea, while for r « T the two quantities differ by a numerical coefficient of 
order unity. In Fig. both I max /e and a are plotted together, and are seen to differ by less 
than a factor of two over the whole T, T range. 

The trapping effect does not happen if the bias is inverted, so that the drain reservoir 
becomes the source and vice versa. In that case we find for the steady-state current the 
expression 

j = 4err 2 (2r + r ) 



r(r + r»(r + 2r ) + 4T 2 (6r + 5r ) 
8err 2 (r 2 + 24T 2 )- 1 if r < r, t 



2eT 2 /r ifr »r,r. (20) 

For strong decoherence the current is the same in both bias directions, but for weak decoher- 
ence the current in the case of inverted bias does not drop to zero but saturates at a finite 
value. The two cases are compared in Fig. [21 We see that the appearance of a maximum 
current as a function of is characteristic for the trapping effect. 

We have for simplicity assumed that all three dots have the same tunnel rates and deco- 
herence rates, but this assumption may be easily relaxed. Let us consider first the case that 
the three-dot structure still has a reflection symmetry, so that dots A and B are equivalent, 
but that dot C has a different tunnel rate V into the reservoir and a different decoherence 
rate T'^. We denote T$ — {T^ + T 1 ^) /2. The result ltTTj|l for the steady-state current generalizes 
to 

4er'T 2 

6 (r*/i> + 2T 2 /r2) 

(21) 



2 + 2T 2 (6 + v/T) + 2r'r (r /r + 2T 2 /r 2 

eT^ if r — > 0, 

2eT 2 /r if r -> oo. 



The steady-state current still contains the desired information on the rates of decoherence, 
with the regimes of weak and strong decoherence governed by Ta, and f^, respectively. 
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Fig. 3 - Solid curves: Dependence of the steady-state current I on the decoherence rate for three 
values of T/T, calculated from Eq. 1181 . The height of the maximum depends non-monotonically on 
F, first increasing oc F and then decreasing oc 1/r, according to Eq. 1191 . Dashed curve: Steady-state 
current if source and drain reservoir are interchanged, calculated from Eq. 121 It for T/T — 20. 



In the most general case of arbitrarily different tunnel rates Ta, Tb, Ta, ~T~b, Tc and deco- 
herence rates T^ A , T^, B , T ( f >c , the steady state current in the limit of weak and strong deco- 
herence takes the form 

i -> ^«( r »x + r» B ) . f r ^ ^ 0j (22a) 



WA + W B 

i — * ; s II 1 d 

waT^ a + WbT^^b + (WA + WBjT^ 



(22b) 



with weight factors 



TaTb FaTb/Ta T b Ta/T b 

wa = 2f+3f ' WA = TXTTV WB = (23) 



that are functions of the tunnel rates — but independent of the decoherence rates. Notice 
that in this asymmetric case the trapped state W wqTb / Ta | A) — ^JwqTa/Tb\B) has unequal 
weights on the two dots A and B. 

In conclusion, we have demonstrated how the well known concept of coherent population 
trapping in atoms may be transferred to a purely electronic system. A large voltage bias plays 
the role of the laser illumination and single-electron tunneling between quantum dots plays the 
role of intra-atomic transitions. Because the quantum dots are charged, the trapped electronic 
state is sensitive to decoherence by coupling to charges in the environment. This decoherence 
destabilizes the trapped state, causing a leakage current / to flow through the quantum dots. 
We have found that the maximal / in the presence of decoherence is set by the trapping rate a, 
with imax ~ ea within a factor of two over the whole parameter range. For small decoherence 
rate Tj, we find I = eT^, which provides a way to measure the coherence time of a charge 
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qubit in a transport experiment. We finally note that extensions to many-electron trapping 
can serve as a source for the formation of entangled electron pairs [20] . 

* * * 

We have benefited from discussions with B. Trauzettel. This work was supported by the 
Dutch Science Foundation NWO/FOM and the US NSF project DMR 0403465. 
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